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Hypcrvelocity  impacts  by  space  particles,  such  as  meteoroids  and  debris,  pose  hazards  to  spacecraft.  The  limits 
of  velocity  of  meteoroid  and  debris  are  derived.  Characteristic  properties  of  hypervelocity  impacts  are  momentum 
transfer,  penetration,  plasma  production,  localization,  and  suddenness.  Using  McDonnell’s  empirical  formulas 
derived  from  laboratory  experiments,  impact  penetrations  and  plasma  production  rates  in  the  space  environment 
arc  calculated.  When  the  critical  temperature  theorem  for  Maxwellian  space  plasmas  is  used,  the  energy  of  the 
plasma  generated  is  shown  to  be  too  low  to  induce  any  significant  spacecraft  charging.  The  plasma  generated, 
however,  can  induce  a  transient,  sustained  or  avalanche  discharge  between  differentially  charged  surfaces.  The 
discharge  current  depends  not  only  on  the  plasma  density  generated  but  also  on  the  neutral  gas  released  on  impact. 
A  scenario  of  impact  induced  hazard  following  days  of  passage  of  a  high-energy  plasma  cloud,  such  as  a  coronal 
mass  ejection  cloud,  is  discussed.  Some  mitigation  methods  are  discussed.  Finally,  we  discuss  whether  electrons 
can  be  accelerated  to  high  energies  in  a  meteor  trail. 
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Nomenclature 

area, 

projectile  diameter,  m 
normal  component  of  meteor  flux,  m^  •  s”' 
electron  distribution  as  a  function  of  electron  energy  £ 
moment  of  inertia,  g  •  m^ 

brightness  in  the  astronomical  scale  (logarithm  of  the 

luminosity  flux  from  a  source) 

spacecraft  mass,  g 

meteor  mass,  g 

electron  mass,  g 

probability  of  meteor  impact 

probability  of  N  impacts 

penetration  depth,  m 

charge  production  rate 

spacecraft  radius,  m 

electron  temperature,  keV 

critical  or  threshold  temperature  for  the  onset  of 

spacecraft  charging,  keV 

time,  s 

projectile  velocity,  km  •  s"' 
relative  velocity,  km/s 
orbital  velocity,  km/s 
escape  velocity,  km/s 
meteor  velocity,  km/s 

yield  of  electrons  emitted  by  a  hypcrvelocity 
particle  impact 

yield  of  ions  emitted  by  a  hypcrvelocity  particle  impact 

number  of  ionizations  generated  by  an  electron 

traveling  through  a  distance  drr 

number  of  electrons  generated  by  an  ion  impact 

on  a  cathode 
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coefficient  of  secondary-electron  emission  with 

primary  electron  energy  £ 

coefficient  of  backscattered-clectron  emission  with 

primary  electron  energy  £ 

number  of  ionizations  generated  by  an  electron 

traveling  from  the  cathode  to  the  anode 

projectile  density,  g  ■  cm“^ 

target  density,  g  •  cm“^ 

electron  impact  ionization  cross  section,  cm”^ 
tensile  strength  of  aluminum 
tensile  strength  of  target 
angular  velocity,  rad/s 


Introduction 

For  any  satellite,  solar  radiation  and  emissions,  space  plasma, 
energetic  particles,  and  meteoroids  provide  sources  of  con¬ 
cern.  Solar  effects,*  space  plasma,^  and  energetic  charged  particle 
interactions^  have  been  addressed  in  many  papers  over  a  period 
of  years.  Hypervelocity  impacts  by  meteoroids  and  space  debris 
on  satellites  have  also  been  studied  but  mostly  for  astronomical 
applications.'' 

The  latter  topic,  that  is,  hypcrvelocity  impact,  is  the  subject  of 
this  paper.  Meteor  impact  may  be  separated  into  two  distinct  parts: 
1)  the  probability  of  collision  between  a  meteoroid  and  a  satellite 
and  2)  the  effects  arising  from  a  collision  between  a  meteoroid 
and  a  satellite.  The  probability  that  a  meteoroid  will  collide  with  a 
satellite  is  proportional  to  the  geometric  cross  section  of  the  satellite 
and  to  the  flux  of  meteoroids  and  other  cosmic  dust  particles.  A 
number  of  authors  have  calculated  such  probabilities^"®  by  taking 
into  account  meteoroid  and  dust  fluxes.**  Three  types  of  interactions 
will  be  considered  in  the  direct  collisions  with  satellite  surfaces:  1) 
penetration  of  the  satellite  surfaces,  2)  collision  with  the  surfaces, 
with  consequent  generation  of  a  cloud  of  neutral  atoms  and  dust, 
and  3)  the  formation  of  plasma  on  collision  of  the  meteoroid  with 
the  satellite  surface  with  the  potential  consequence  of  a  sudden 
discharge. 

The  long-duration  exposure  facility  (LDEF)  had  a  six -year  space 
sojourn,  resulting  in  a  treasure  trove  of  debris-  and  meteoroid-impact 
data.  Studies  of  craters  on  a  variety  of  surfaces  made  possible  a  deter¬ 
mination  of  the  flux  of  particles  as  well  as  a  qualitative  discrimina¬ 
tion  between  space  debris  and  interplanetary  dust  or  meteoroids.  ' 
The  penetration  depth  of  hypervelocity  impacts  is  found  to  have 
an  empirical  relationship  that  depends  on  the  relative  densities  of 
the  meteoroids  and  target,  the  meteoroid  diameter,  and  the  impact 
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Table  1  Spacecraft  failed  and  damaged  due  to  meteoroid  or  debris  impact^ 


Vehicle 

Date 

Effect 

Possible  cause 

HST 

April  1990 

5000  impacts  in  4  years, 
solar  cell  punctured. 

Meteoroid/debris 

ISEE-1 

Oct.  1977 

Detector  window  punctured. 

Meteoroid 

lSEE-1 

Aug.  1978 

All  isobutane  gas  lost  in  5  days. 

Meteoroid 

Olympus 

Aug.  1993 

Satellite  failed  to  function. 

Meteoroid 

SEDS  2 

March  1994 

Mission  terminated. 

Meteoroid/debris 

MSTI  2 

March  1994 

Satellite  communication  lost. 

Meteoroid/debris 

Kosmos-1275 

July  1981 

Fragmented  into  200  pieces. 

Debris 

Solar-A 

Aug.  1991 

Telescope  punctured. 

Meteoroid 

STS-45 

March  1992 

Gouges  on  wing  edge. 

Debris 

STS-49 

May  1992 

Chip  in  window  pane. 

Debris 

“Data  adapted  from  Koons  et  al.^*  HST,  Hubble  Space  Telccopc;  ISEE,  International  Sun-Earth  Explorer;  and 
STS,  Space  Transportation  System. 


velocity.'^' Laboratory  simulations  of  the  hypervclocity  impact  of 
debris  on  surfaces  have  been  performed  extensively.  For  example, 
in  the  impact  of  a  70-nm  boron  carbide  particle  (4  x  10“’^  g)  at  a 
speed  of  94  km  •  s“'  on  an  aluminum  target  doped  with  silver,  the 
total  energy  of  the  products  (vaporization,  ionization  and  kinetic  en¬ 
ergy  of  products)  represented  about  3.7%  of  the  projectile  energy.''’ 
In  addition  to  meteoroids,  space  debris  represent  another  type  of 
interaction  that  somewhat  resembles  the  behavior  of  meteoroids. 
Because  space  debris  is  for  the  most  part  banded  at  an  altitude 
near  1000  km,  this  interaction  can  be  considered  to  be  important 
in  the  altitude  region  600-1600  km  (Ref.  15).  This  hazard  arises 
from  the  breakup  of  satellite  and  spent  upper  stages,  leading  to  the 
presence  of  large  amounts  of  debris.  It  is  estimated  that  the  number 
of  space  debris  objects  greater  than  0.01  m  diam  will  increase  from 
near  8000  in  2000  to  ~180,000  in  2100  at  geosynchronous  orbit.’^ 
In  view  of  the  rapidly  increasing  number  of  spacecraft  and  the  large 
increase  in  the  number  of  space  debris  objects  in  the  next  century, 
hypervelocity  impact  on  spacecraft  systems  needs  to  be  considered 
as  an  important  component  of  spacecraft  interactions  and  spacecraft 
survivability. 

Hypervelocity  Impacts  as  Space  Hazards 

When  meteoroids  enter  the  atmosphere,  they  do  so  at  velocities 
ranging  from  about  11  to  72  km  - s"’.  Space  debris  have  orbital 
velocities  of  7-1 1  km  •  s“'  typically.  These  large  velocities  imply 
large  kinetic  energies  and  momenta. 

For  known  events,  we  refer  to  the  Koons  et  al.'*  list  of  missions 
lost  or  terminated  due  to  the  space  environment  from  1973  to  1997. 
Of  the  1 0  cases  listed  in  Table  1 ,  four  were  due  to  meteoroid  impacts, 
whereas  the  rest  were  attributed  to  debris  or  meteoroids/debris.  The 
three  spacecraft  lost  due  to  meteoroid  impact  are  Olympus  (August 

1993) ,  the  Small  Expendable  Deployer  System  (SEDS)  (2  March 

1994) ,  and  the  Miniature  Sensor  Technology  Integration  (MSTI)  (2 
March  1994).'®  Some  of  the  electrostatic  discharge  cases  could  be 
due  to  meteor  impact  as  well.  Therefore,  meteor  and  debris  impacts 
on  spacecraft  constitute  a  considerable  percentage  of  the  total  cases 
of  space  missions  lost. 

Whereas  spacecraft  interactions  with  space  plasmas  have  been 
well  studied,^’ our  understanding  of  hypervelocity  impacts  on 
spacecraft  arc  in  a  less  mature  stage  of  development.  In  this  pa¬ 
per,  meteoroid  velocities,  impact  penetration  depth,  impact  plasma 
production,  and  hazards  of  hypervelocity  impacts  are  studied. 

The  organization  of  this  paper  is  as  follows.  The  velocities  of  me¬ 
teors  and  the  kinetic  energies  of  atoms  at  such  velocities  are  brack¬ 
eted.  Because  of  the  tremendous  impact  kinetic  energies,  significant 
surface  penetration  may  occur.  Impact  penetration  depths  are  cal¬ 
culated  using  empirical  formulas  based  on  experiment.  Two  likely 
effects  due  to  hypervelocity  impact,  namely,  mechanical  and  elec¬ 
trical,  are  studied.  Mechanical  effects  include  surface  or  structural 
damage  by  impact  penetration  and  perturbation  of  spacecraft  linear 
and  angular  momenta.  Electrical  effects  are  due  to  plasma  generation 
on  impact.  A  most  common  concern  is  whether  the  impact-generated 
plasma  can  cause  spacecraft  charging  and  discharging.  After  exam¬ 
ining  the  laboratory  results  for  the  energy  of  the  impact-generated 
plasma,  we  conclude  that  meteoric  impact  is  unlikely  to  cause  sig¬ 
nificant  spacecraft  charging.  The  plasma,  however,  disturbs  scien- 
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Fig,  1  Orbital  and  escape  velocities  at  various  altitudes. 

tific  measurements.  More  significantly,  the  plasma  can  cause  current 
leakage  or  even  avalanche  discharges  between  differentially  charged 
configurations.  Some  most  likely  differentially  charged  configura¬ 
tions  are  solar  cells,  surfaces  of  very  different  properties,  and  double¬ 
layer  deep  dielectric  charging.^"  These  unstable  configurations  pose 
significant  potential  hazard  for  spacecraft.  Another  hazard  involving 
both  mechanical  and  electrical  aspects  is  pointed  out.  If  an  impact 
penetrates  through  a  wall,  the  impact-produced  plasma  can  exit  from 
both  sides  of  the  puncture.  The  plasma  inside  the  wall  may  cause 
significant  short  circuits  among  the  electronic  components  inside, 
this  being  a  significant  hazardous  scenario.  Finally,  some  mitigation 
techniques  are  discussed. 

Meteoroid  and  Debris  Velocities 

The  escape  velocity  Vs  of  a  particle  from  a  point  in  the  solar 
system  and  the  orbital  velocity  V„  at  the  point  are  related  by 

Vl  =  2V^  (1) 

A  meteoroid  at  the  Earth’s  location  has  an  escape  velocity 
V(s  =  42.0  km/s.  For  its  head-on  collision  with  the  Earth,  which 
has  orbital  velocity  of  29.7  km/s,  the  relative  velocity  is  42.0-1- 
29.7  =  71.7  km/s.  Including  the  Earth’s  gravitation,  the  maximum 
relative  velocity  is  72.6  km/s. 

The  velocity  needed  to  escape  from  Earth  from  an  altitude  of 
120  km  is  about  11.1  km/s.  Because  the  orbital  velocity  is  related 
to  the  escape  velocity  [Eq.  (1)]  the  orbital  velocity  of  debris  at 
120  km  is  ll.l/.y2  =  7.8  Im/s  approximately  (Fig.  1).  Below  that 
altitude,  atmospheric  drag  becomes  significant  in  the  deceleration 
of  meteorites  and  particulates. 

If  the  eccentricity  of  the  Earth’s  orbit  is  included,  the  Earth’s 
orbital  velocity  ranges  from  29.3  to  30.3  km/s,  with  heliocentric 
escape  velocities  of  41.4-42.8  km/s.  As  a  result,  the  impact  speeds 
could  range  up  to  73.1  km/s  or,  when  the  gravitational  attraction  is 
taken  into  account,  up  to  74.0  km/s. 

The  preceding  treatment  assumes  that  all  meteoroids  are  on  closed 
solar  system  orbits.  There  is  now  conclusive  evidence  for  a  signifi¬ 
cant  population  of  meteoroids  that  arrive  from  interstellar  space.  At 
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SPACE  ENVIRONMENT  EFFECTS 


Impact  Penetration,  Puncture, 

Debris  &  Meteorites 

Plasma  Production 

Single  Event  Upsets 

Cosmic  Rays 

GeV 

Solar  Eruptions 

Shocks,  Geomagnetic  Storms 

Radiation  Belt  Particles 

Radiation  Damage,  Degradation, 

MeV 

Deep  Dielectric  Charging 

Energetic  Plasma 

Surface  Charging  &  Discharging, 

Spacecraft  Anomalies 

keV 

Low*energy  Plasma 

Leakage,  Erosion 

eV 

Neutral  0  Atoms 

Chemical  Reactions,  Erosion 

Fig.  2  Energy  scale  of  various  spacecraft-environment  interactions. 


the  mass  range  10  kg,  they  probably  account  for  a  few  percent  of 
the  total  meteoroid  input.  There  are  some  radar  records^'  indicating 
a  tail  of  higher  velocity  meteors.  Relatively  recent  reviews  of  the 
topic  of  detection  of  meteoroids  with  origin  outside  the  solar  system 
arc  given  by  Hawkes  ct  al.^^  and  by  Ceplecha  ct  al.^^ 

Impact  Energy 

Because  of  the  largc-impacl-energy,  hypcrvciocity  impacts  of  me¬ 
teoroids  pose  hazards  to  space  systems^'*  and  human  space  flights.^^ 
A  Leonid  meteoroid  of  0. 1  g  in  mass,  traveling  at  7 1  km/s,  would  im¬ 
pact  the  solar  cell  of  a  spacecraft  with  a  force  equal  to  that  of  a  bullet. 
A  calcium  atom,  for  example,  traveling  at  72  km/s  has  a  kinetic  en¬ 
ergy  of  about  1  kcV.  Therefore,  a  meteoric  particle,  often  with  more 
than  millions  of  atoms,  has  a  kinetic  energy  exceeding  1  GeV.  Both 
debris  and  meteoric  particles  can  cause  physieal  damage.  Debris 
and  mcleorie  impacts  arc  the  most  energetic  interactions  compared 
with  all  others,  such  as  impacts  by  neutral  oxygen  atoms  in  the  iono¬ 
sphere,  hot  plasmas  at  geosynchronous  orbits,  energetic  electrons 
and  ions  in  the  radiation  belts,  coronal  mass  ejection  particles,  and 
cosmic  rays  (Fig.  2). 

Fast  atoms  and  ions  can  penetrate  into  surfaces  on  impact  gener¬ 
ating  ionization  tracks  inside  the  solid,  ejecting  plasma  and  neutral 
gas  vapor  to  outside  the  surface,  and,  if  puncture  occurs,  spraying 
plasma  and  energetic  projectiles  to  behind  the  solid  wall.  Some  of 
the  hazardous  effects  are  discussed  in  a  later  section. 

Impact  Penetration 

Cour-Palais^^  measured  in  the  laboratory  the  penetration  depth  p 
as  a  function  of  the  impacting  particle  diameter  d,  density  p,  and 
velocity  V  and  obtained  an  empirical  equation, 

=  (2) 

where  p„  is  in  grams  per  cubic  centimeter,  V  is  in  kilometers  per 
second,  and  C  is  a  constant  that  is  characteristic  of  the  target  ma¬ 
terial  and  its  condition.  An  approximate  formula  based  on  LDEF 
observations  is  given  in  Tribble^’: 

p  =  (3) 

where  m  is  the  mass  (grams)  of  the  panicle,  pr  the  target  density 
(grams  per  cubic  centimeter),  V  the  normal  component  of  the  par¬ 
ticle  velocity  relative  to  the  surface,  and  K  a  materials  constant. 

McDonnell  and  Sullivan"  and  McDonnell  et  al.'^  obtained  an 
empirical  equation  for  meteoroids  and  debris  impacts  applicable 
under  a  wide  range  of  target  densities  pr  and  tensile  strengths  ot- 

p/d  =  0.7658cl‘’“'^(Pp/pr)‘’‘"'*((TAi/oT)‘’’^''v‘’®“  (4) 

Note  that  these  formulas  are  approximate  because  they  are  based  on 
laboratory  experiments  with  solid  particles.  In  space,  the  meteoroids 
arc  likely  to  be  fluffy  and  irregularly  shaped^®;  hence,  the  actual 


densities  and  tensile  strengths  may  be  incorrectly  represented  in  the 
preceding  formulas. 

Using  the  McDonnell  et  al.'^  formula  [Eq.  (4)],  we  have  calcu¬ 
lated  the  penetration  depth  to  be  0.35  mm  for  a  1.3  x  10“^  g  Leonid 
meteoroid  colliding  with  a  typical  satellite  mylar  surface  at  71  km/s. 
A  penetration  of  0.35  mm  into  typical  solar  panel  Mylar®  surface 
would  cause  damage  on  spacecraft  solar  panels.  The  effect  of  such 
a  penetration  may  put  out  a  single  solar  cell,  leading  to  degrada¬ 
tion  of  performance  over  time.  A  2.9  x  10  '*  g  Leonid  meteoroid 
would  penetrate  0.01  m  into  aluminum.  It  is  difficult  to  predict  how 
damaging  an  impact  penetration  would  be  because  it  would  depend 
on  the  details  of  what  systems  were  near  the  impact.  Penetrations 
followed  by  plasma  and  current  generations,  to  be  discussed  later, 
would  be  more  hazardous,  affecting  wider  areas. 

A  comparison  of  formulas  (3)  and  (4)  is  now  discussed.  The 
penetration  depth  p  obtained  by  using  the  LDEF  formula  [Eq.  (3)]  is 
0.09  m  for  the  2.9(— 7)  kg  particle.  This  is  compared  to  0.01  m  using 
McDonnell’s  formula  [Eq.  (4)].  The  problem  with  using  LDEF  data 
is  that  they  do  not  include  any  impacts  from  particles  this  large.  Note 
the  small  probability  of  such  an  impact  even  in  a  very  large  Leonid 
storm.  LDEF  really  stops  at  about  1  (—8)  kg  as  far  as  anything  that 
is  statistically  meaningful. 

Impact  Probability 

The  probability  P  of  impact  by  a  particle  on  a  surface  is  given  by 

/(  + A/ 

ditF(t)A  (5) 

where  At  is  the  duration  of  the  exposure  and  F{t)  the  normal  com¬ 
ponent  of  the  flux  of  particles.  To  calculate  the  probability  of  a 
given  impact  penetration  depth  p,  one  uses  McDonnell’s  formula 
[Eq.  (4)]  to  get  a  minium  diameter  d  to  cause  the  event.  With  an 
assumed  particle  density,  this  gives  the  mass  m.  An  empirical  re¬ 
lationship  between  mass  m  and  brightness  M  on  the  astronomical 
scale,  which  is  a  logarithmic  scale  of  luminous  flux  from  a  source, 
is  given  by^^ 

logio  m  =  -0.4M  -  1 0.97  -F  1 .7  log,o  V  (6) 

where  V  is  in  centimeters  per  second.  Then  the  meteor  shower 
distribution  function,  a  function  of  brightness  M  or  mass  m,  gives 
one  the  flux  Fit).  That,  along  with  the  area  A,  gives  the  probability 
[Eq.  (5)]. 

The  slope  of  the  meteor  shower  distribution  function  is  commonly 
characterized  by  the  so-called  population  index  r,  which  gives  the 
relative  number  of  meteors  detected  in  successive  divisions  of  unit 
magnitude  M .  Equation  (6)  is  applied  to  meteor  mass  at  the  limiting 
magnitude  6.5M.  Details  on  modeling  Leonids  meteor  shower  pop¬ 
ulation  distribution  as  a  function  of  meteor  mass  is  given  in  Ref.  30. 

Whereas  Eq.  (6)  has  been  used  for  over  40  years,  more  recent 
observational  work^"^^  seems  to  suggest  slightly  different  values  of 
the  numerical  factors.  Recent  studies^^"^®  have  suggested  different 
mass-magnitude  relationships.  For  our  purposes,  however,  Eq.  (6) 
seems  adequate. 

The  meteoroid  flux  is  expected  to  be  unchanged  from  the  entry 
point  down  to  about  120  km,  where  meteoric  ablation  due  to  fric¬ 
tional  heating  begins.  We  examine  two  cases,  the  penetration  of  an 
unshielded  solar  cell  and  the  penetration  of  a  heavily  shielded  space¬ 
craft  surface.  For  the  first  case,  we  model  the  penetration  as  a  5-mm 
pit  in  a  Mylar  surface.  The  smallest  Leonid  needed  to  form  this  pit 
is  0.35  mm  in  diam  and  weighs  1.3  x  10“^  g.  For  a  major  Leonid 
storm  such  as  that  of  1966,  the  flux  of  such  particles  is  8.63  x  10  ® 
m“^s“'  and,  assuming  20  m^  of  surface  area  and  a  2.5-h  storm,  the 
probability  of  such  a  hit  is  1.3%  according  to  Eq.  (5).  In  the  sec¬ 
ond  case,  we  assume  a  0.05-m  pit  in  an  aluminum  surface,  which 
requires  a  Leonid  meteor  of  4.4-mm  diam  or  0.028  g.  Again  for  a 
major  Leonid  storm,  the  flux  would  be  4.69  x  10^*°  m^^  s"'  and 
the  probability  is  0.008%. 

From  elementary  probability  theory,  the  probability  Pn  of  N  im¬ 
pacts,  each  one  of  which  is  uninfluenced  by  the  other,  in  a  duration 
At  is  given  by  the  Poisson  distribution: 

Pn  =  P"exp(-P)/A! 


(7) 
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If  W  =  1  and  P  <K  1 ,  one  recovers  Eq.  (5)  from  Eq.  (7)  by  using  the 
Taylor  expansion  of  the  exponential. 

Impact  Probabilities  of  Annual  Meteor  Showers 

In  the  last  section,  some  probabilities  were  quoted  for  penetra¬ 
tion  of  a  solar  panel  surface  and  a  heavily  shielded  0.05-m  alu¬ 
minum  panel  during  a  Leonid  storm  of  historic  magnitude.  Here, 
these  probabilities  arc  compared  to  those  encountered  routinely  by 
spacecraft  during  the  annual  meteor  showers. Three  impact  pene¬ 
tration  scenarios  are  defined,  and  for  each  material,  a  critical  pene¬ 
tration  depth  and  a  surface  area  exposed  to  the  meteor  stream  need 
to  be  defined.  First,  a  5-mm  penetration  into  a  Mylar  surface,  which 
would  most  likely  shatter  a  solar  panel  segment,  is  chosen.  Using  the 
DSCS  satellite  as  typical,  the  exposed  surface  area  of  solar  panels 
is  around  20  m^.  The  second  case  is  for  a  lightly  shielded  surface  of 
the  spacecraft.  This  could  be,  for  example,  an  antenna  section  or  a 
surface  containing  scientific  instruments.  Here  a  critical  penetration 
of  10  mm  in  aluminum  and  a  surface  area  of  10  m^  are  assumed. 
Only  one  face  of  the  satellite  would  be  exposed  to  the  flux  at  any 
given  time.  For  the  third  case,  the  equivalent  of  the  critical  space¬ 
craft  design  shield  (USDS),®  which  is  defined  by  a  49-mm  pit  in 
aluminum,  is  selected.  The  exposed  area  is  taken  to  be  10  m^  as 
typical  of  a  medium-sized  spacecraft. 

We  show  in  Fig.  3  the  impact  rates  of  the  annual  meteor  show¬ 
ers  for  an  equivalent  penetration  of  one  CSDS  unit  in  aluminum, 
throughout  the  year,  and,  similarly,  in  Fig.  4,  the  impact  rates  for 
a  1-mm  penetration  in  aluminum.  In  terms  of  maximum  rate,  the 
storm  Leonids  (LEO  in  Fig.  4)  are  the  most  likely,  but  the  duration 
of  the  Leonids  is  relatively  short  compared  to,  for  example,  the  rj- 
Aquarids  (ETA  in  Fig.  4).  It  is  the  combination  of  peak  intensity 
and  duration  that  ultimately  determines  the  probability  of  an  impact 


Fig.  3  Annual  distribution  of  the  probability  per  square  meter  per 
second  of  meteoroid  impact  to  1  CSDS  unit  (4.9  mm)  in  aluminum. 


Fig.  4  Annual  distribution  of  the  probability  per  square  meter  per 
second  of  meteoroid  impact  to  1  mm  in  aluminum. 


Table  2  Some  major  showers  and  impact  probabilities 


Shower 

Characteristics 

Impact  probability,  % 

Dates 

ZHR 

V 

Case  1 

Case  2 

Case  3 

Leonid  storm 

17  Nov. 

150,000 

71 

2.00 

0.067 

0.0036 

Annual  leonids 

14  Nov.-21  Nov. 

23 

71 

0.48 

0.0086 

0.00003 

i)-Aquarids 

1 9  April-28  May 

37 

66 

2.32 

0.035 

0.0002 

Perseids 

17  July-24  Aug. 

84 

61 

4.39 

0.058 

0.0003 

Orinids 

2  Oct.-7  Nov. 

25 

67 

2.48 

0.029 

0.0001 

from  a  particular  shower.  We  see  that  the  probability  for  impact  is 
distributed  pretty  much  uniformly  throughout  the  year,  with  a  rela¬ 
tive  lull  in  January-March  (solar  longitude  >270°).  Parameters  for 
the  showers  are  taken  from  Ref.  37. 

In  Table  2  we  have  compared  probabilities  for  the  three  model 
impact  scenarios  for  a  major  Leonid  storm  and  for  four  of  the  most 
intense  annual  meteor  showers.  The  mass  distributions  used  to  com¬ 
pute  the  flux  of  the  limiting  mass  required  to  inflict  the  penetration 
were  computed  as  described  by  McNeil  et  al.^°  The  ranges  of  dates 
chosen  for  the  calculation  are  those  listed  for  these  showers  by 
the  International  Meteor  Organization’s  Meteor  Shower  Calendar 
(available  on  line  at  URL:  www.imo.net).  Note  that,  in  computing 
the  probabilities  in  Table  2,  a  relatively  complete  list  of  50  major 
annual  showers,  compiled  by  Jenniskens,^’  was  used,  and  therefore, 
the  total  probabilities  contain  contributions  from  other  less  intense 
showers  during  these  time  periods.  Even  so,  this  represents  the  total 
flux  that  a  satellite  would  encounter  during  the  time  periods  listed. 

Note  that,  for  cases  1  and  2,  penetration  of  0.005  m  of  Mylar 
and  10  mm  of  aluminum,  the  probabilities  from  the  annual  showers 
are  comparable  to  and  often  greater  than  those  from  a  truly  historic 
meteor  storm  like  the  Leonids  in  1966.  This  result  arises  from  two 
factors.  First,  the  storm  is  of  short  duration,  a  matter  of  hours,  while 
some  annual  showers  last  from  a  few  weeks  to  a  month,  albeit  not 
at  their  maximum  zenith  hourly  rate  (ZHR).  Second,  the  mass  dis¬ 
tributions  of  the  annual  showers  favor  smaller  particles,  whereas 
the  storm,  with  its  relatively  fresh  cometary  ejecta,^®  has  a  dispro¬ 
portionately  large  number  of  larger  meteors.  It  is  only  for  the  truly 
catastrophic  impacts,  such  as  case  3,  where  the  Leonid  storm  poses 
more  of  a  hazard,  and  even  then,  only  with  a  very  small  probability 
and  only  a  factor  of  10  greater  than  the  annual  showers. 

The  probabilities  for  penetrating  a  solar  panel,  case  1,  are  rela¬ 
tively  large  and,  over  the  course  of  a  year,  add  up  to  one  chance  in 
five  or  so.  Over  a  year,  the  probability  of  penetration  of  a  lightly 
shielded  surface  on  a  single  satellite  is  about  1  in  1000.  This  may 
seem  small,  but  assuming  1000  satellites  in  orbit,  such  a  penetration 
on  at  least  one  of  them  becomes  likely.  The  probability  of  pene¬ 
tration  of  a  heavily  shielded  surface,  even  during  a  Leonid  storm, 
remains  small.  However,  note  that  this  increases  linearly  with  ex¬ 
posed  area.  For  a  very  large  object,  such  as  a  space  station,  with 
1000-m^  exposed  area,  the  probability  of  a  0.0 1-m  penetration  rises 
to  6%  during  a  Leonid  storm,  with  the  probability  for  a  CSDS  pen¬ 
etration  at  0.36%.  These  are  not  negligible,  especially  on  crewed 
spacecraft.  Finally,  note  that,  in  general,  satellites  are  not  designed 
with  the  CSDS  level  of  protection,  so  that  damage  may  certainly 
take  place  with  impacts  from  particles  significantly  smaller  than 
those  implied  by  this  limit. 

Perturbation  of  Angular  Momentum 

To  estimate  the  effect  on  angular  momentum,  a  spherical  space¬ 
craft  of  mass  A/,  =  10®  g, radius/?  =  2  m,  and  a  solar  panel  extending 
10  m  is  considered.  A  Leonid  meteor  of  velocity  u  =72  km/s  and 
mass  m  =  0. 1  g,  colliding  with  the  tip  of  a  boom  or  an  extended 
solar  panel  would  add  an  angular  velocity  Aw  given  by 

I  Aw  =  mvr  (8) 

where  I  is  the  moment  of  inertia  of  the  spacecraft,  m  =  0.1  g, 
v  =  12  km/s,  and  r  =  12  m.  For  a  sphere,  the  moment  of  iner¬ 
tia  is  given  by  /  =  The  result  is  Aw  =  0.54  rad  •  s^’  or 

0.086  Hz.  This  magnitude  of  spin  rate  change,  by  itself,  is  not  sig¬ 
nificant  and  would  not  disrupt  spacecraft  operations.  Nevertheless, 
it  is  detectable  in  certain  telemetry  signals  such  as  electric  field  mea¬ 
surements  by  means  of  booms  rotating  with  the  spacecraft.  One  way 


no 
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to  differentiate  a  meteor  impact  from  one  by  a  debris  is  to  look  at  the 
vector  direction  of  the  sudden  change  of  spacecraft  angular  momen¬ 
tum,  from  which  one  deduces  the  direction  of  impact.  Typically,  de¬ 
bris  have  near  circular  orbits,  depending  on  the  satellite  from  which 
the  debris  comes,  and  arc  much  heavier  but  slower  (at  15  km  •  s^' 
for  head-on  collisions  at  low  Earth  orbits)  than  meteoroids. 

Total  Charge  Production 

When  a  meteor  collides  with  a  surface,  the  collision  energy  is 
enough  to  generate  a  plasma  cloud  on  the  surface.  McDonnell 
cl  al.'^  obtained  from  extensive  laboratory  experiments  an  empirical 
formula  for  the  electron  charge  production  rate  Q ,  when  a  meteoroid 
collides  with  an  aluminum  surface: 

Q  =  10-°®%‘°^(l//5)^''*  (9) 

The  charge  production  rate  is  calculated  for  minor,  moderate,  and 
major  Leonid  showers  of  ZHR  =  100,  5000,  and  150,000,  respec¬ 
tively.  Assuming  a  population  index  of  1.8  for  the  Leonids^®  one 
obtains  the  time-averaged  charge  production  rales  during  a  3-h 
storm  period  as  0.2  x  10  1.0  x  10^^,  and  3.0  x  10“^  nA/cm'^. 

For  comparison,  the  average  quiet  lime  ambient  flux  of  charged 
particles  on  SCATHA  at  near  geosynchronous  altitudes  was  calcu¬ 
lated  by  Purvis  et  al.‘'°  to  be  0.1 15  nA/cm^.  Similar  results  of  quiet 
time  geosynchronous  altitude  ambient  flux  have  been  obtained  by 
using  the  current-voltage  cutoff  point  of  electron  beam  emission 
from  SCATHA.'”  Thus,  the  time-averaged  meteor  impact  generated 
charge  flux  is  less  than  the  ambient  flux  at  quiet  times  on  SCATHA. 
According  to  this  analysis,  the  time-averaged  meteor  flux  would  not 
induce  significant  charging.  However,  meteor  impacts  arc  different 
from  average  flux  both  because  they  arc  localized  and  nearly  in¬ 
stantaneous.  It  is  difficult  to  say  whether  these  differenees  would 
cause  unwanted  effects.  There  are  other  independent  reasons  why 
charging  by  meteor  impact  is  unlikely,  however,  which  we  discuss 
in  the  next  section. 

Can  Meteor  Impact  Induce  Spacecraft  Charging? 

It  is  common  to  associate  the  rate  of  plasma  production  by  meteor 
impact  with  the  likelihood  of  spacecraft  charging,  although  that 
association  has  never  been  established.  If  meteor  impact  can  induce 
hundreds  or  thousands  of  volts  on  spacecraft,  it  would  be  significant. 
We  show  next  that  this  is  very  unlikely. 

Before  we  discuss  theory,  we  show  an  interesting  observation 
(Fig.  5)  that  tempts  one  to  associate  instantaneous  spacecraft  charg¬ 
ing  with  meteor  impact.  In  Fig.  5,  the  lower  panel  plots  the  kp  index. 
When  kp  is  low  (below  about  4),  the  magnetosphere  is  quiet;  when 
kp  is  high,  the  magnetosphere  is  stormy.  As  Fig.  5  shows,  it  was 


stormy  on  14  November  but  quiet  on  15  and  16  November  1998. 
The  midpanel  plots  the  electron  flux  measured  by  the  GOES  satel¬ 
lite.  The  flux  fluctuates  wildly  during  stormy  periods.  TTie  ambient 
ion  spectrum  (not  shown)  measured  on  DSCS  also  did  not  show 
stormy  or  unusual  behavior  on  16  November.  It  is  well  known  that 
during  severe  geomagnetic  storms  or  substorms,  spacecraft  charging 
sometimes  occurs  if  the  ambient  plasma  becomes  energetic.  There 
was  no  charging  on  the  DSCS  satellite  on  14  and  15  November, 
although  two  brief  instantaneous  events  (about  3000  and  10,000  s 
universal  time)  were  recorded  on  16  November  (top  panel).  Because 
1 6  November  1 988  was  a  quiet  day  (not  stormy),  any  cause  of  charg¬ 
ing  due  to  hot  ambient  plasma  is  ruled  out.  The  peak  of  the  Leonids 
shower  of  1998  arrived  on  16  November  1998  but  not  at  the  time 
of  the  two  charging  events.  It  is  tempting  to  associate  the  charging 
events  with  meteor  impacts,  but  the  true  cause  is  unknown. 

A  basic  theorem  in  spacecraft  charging  states  that  in  a  Maxwellian 
plasma  environment,  there  exists  a  critical,  or  threshold,  plasma  elec¬ 
tron  temperature  for  surface  charging  to  occur.^  ’^  The  reason  for 
this  theorem  is  as  follows.  At  equilibrium,  Kirchhoff’s  circuit  law 
requires  the  surface  potential  0  be  determined  by  a  current  balance 
equation.  In  space  and  in  the  laboratory,  the  electron  current  is  nor¬ 
mally  two  orders'’^  '*^  of  magnitude  greater  than  that  of  ions  because 
ions  are  much  heavier  and  slower.  Because  an  initially  uncharged 
spacecraft  placed  in  a  plasma  would  intercept  more  electrons  than 
ions  because  of  mass  difference,  the  spaceeraft  would  charge  to  a 
negative  potential,  which  renders  the  fluxes  to  balance  at  equilib¬ 
rium. 

Photoelectrons  arc  important  in  sunlight.  Low-level  positive 
charging  occurs  when  the  ambient  plasma  environment  is  quiet 
while  the  surface  is  in  sunlight.'  The  level  is  low  beeause  the  domi¬ 
nant  sunlight  spectral  line,  Lyman-a,  has  about  10.2  eV,  whereas  the 
work  function  of  most  surface  materials  is  about  4-5  eV,  resulting 
in  photoelectron  energies  of  a  few  electron  volts  only.  When  the 
charging  level  exceeds  a  few  electron  volts,  the  photoclectrons  can 
not  leave.  Thus,  low-level  charging  by  sunlight,  at  geosynchronous 
orbits,  for  example,  is  up  to  a  few  electron  volts  only.  However,  if 
the  spacecraft  is  near  the  sun,  the  far  UV  lines  become  more  intense. 
As  a  result,  more  high-energy  (tens  of  electron  volts)  photoelectrons 
will  come  out,  in  contrast  to  the  geosynchronous  orbit  case  where  the 
photoclectrons  are  of  1  to  few  electron  volts  only.  Consequently,  the 
spacecraft  can  charge  to  tens-  of  volts  (positive)  when  approaching 
the  sun,  in  contrast  to  1  or  few  volts  at  the  geosynchronous  orbits. 

We  consider  negative  charging  here  because,  for  the  case  un¬ 
der  discussion,  the  spacecraft  interacts  with  the  plasma  generated 
by  meteoroid  impact.  Charging  by  plasmas  needs  to  be  explained 
clearly.  Electron  flux  is  higher  than  ion  flux  because  of  mass  dif¬ 
ference.  Therefore,  any  uncharged  object  placed  in  a  plasma  would 


Fig.  5  Top:  Two  brief  charging  events  of  kV  magnitude  (negative  voltage)  on  the  DSCS  satellite,  16  Nov.  1998;  none  was  on  14  and  IS  Nov.  Middle: 
GOES  satellite  measurements  showing  moderately  stormy  plasma  condition  on  14  Nov.  but  quiet  condition  on  16  Nov.  Bottom:  geomagnetic  index 
showing  low  ip  on  16  Nov.  (The  GOES  and  kp  data  were  obtained  from  the  Space  Environment  Center,  Boulder,  CO,  National  Oceanic  and  Atmospheric 
Administration,  U.S.  Department  of  Commerce.) 
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intercept  more  electrons  than  ions.  As  a  result,  the  object  becomes 
negatively  charged  to  an  equilibrium  potential,  which  renders  the 
incoming  flux  of  electrons  to  equal  the  ion  flux.  The  rendering  is 
carried  out  by  the  repulsion  (to  electrons)  and  the  attraction  (to  ions) 
of  the  negative  potential. 

To  calculate  the  onset  of  negative  charging,  it  is  often  a  good 
approximation  to  ignore  the  ions  because  the  ion  flux  is  orders  of 
magnitude  smaller  than  the  electron  flux.  The  approximate  current 
balance  is,  therefore,  between  the  incoming  ambient  electrons  and 
the  outgoing  secondary  and  backscattered  electrons: 

poo  poo 

/  dEEf(E)=  /  dEEf(E)lS(E)  +  t](E)]  (10) 

Jo  Jo 

In  Eq.  (10),  the  left-hand  side  represents  the  incoming  electron  flux 
with  an  energy  distribution  /(£);  the  right-hand  side  represents  that 
of  the  outgoing  electrons.  Here,  8(E)  and  t](E)  are  the  secondary 
and  backscattered  electron  coefficients.  It  is  more  convenient  to  use 
energy  integration  than  velocity  integration  because  the  coefficients 
are  measured  as  functions  of  energy.  For  normal  electron  incidence, 
the  angular  variables  cancel  in  Eq.  (10).  For  angle-distributed  in¬ 
coming  flux,  the  coefficients  have  to  be  angle  dependent  and  the 
angular  integrations  have  to  be  included. 

Substituting  the  Maxwellian  distribution  /(£), 

/(£)  =  n(mc/27r^r)5  exp(— £//:£)  (11) 

in  the  current  balance  equation  (10),  one  finds  that  the  density  n 
cancels  on  both  sides.  Thus,  if  the  assumption  is  made  that  the 
energy  distribution  of  the  impact  generated  plasma  is  Maxwellian, 
the  current  balance  [Eq.  (10)]  is  independent  of  the  plasma  density  n 
and  is  a  function  of  electron  temperature  T  only.  This  is  an  important 
theorem.  The  solution  T*  of  Eq.  (10)  is  the  critical,  or  threshold, 
temperature  for  the  onset  of  surface  charging.  Typical  values  of 
T*  are  1000-2000  eV  depending  on  the  surface  material.^’'’’'’^  For 
Si02,  for  example,  T*  is  1.7  keV  for  normal  electron  incidence 
and  2.6  keV  for  isotropic  incidence.  Thus,  if  the  plasma  electron 
temperature  is  below  1.7  keV,  a  silicon  dioxide  surface  charging 
onset  cannot  occur. 

When  70-nm  boron  carbide  projectiles  of  94  km/s  were  used  on 
aluminum  surfaces  in  the  laboratory,  it  was  found  that  2.9%  of  the 
projectile  energy  went  to  the  plasma  kinetic  energy,  0.5%  went  to 
ionizing  the  neutrals  to  form  plasma,  and  0.3%  went  to  vaporizing 
the  neutrals.'*'  The  hydrogen  line  in  the  plasma  showed  120  eV  and 
the  carbon  line  40  eV.  Griin  (personal  communication,  2000)  also 
observed  that  the  impact  of  meteor-sized  projectile  generated  plasma 
energy  of  typically  10  to  a  few  tens  of  electron  volts  only.  In  any  case. 


the  energy  distribution  of  the  impact-generated  plasma  has  not  been 
measured.  If  the  assumption  is  made  that  the  energy  distribution  of 
the  impact-generated  plasma  is  Maxwellian,  the  plasma  temperature 
is  the  key  factor  determining  whether  significant  negative-voltage 
charging  can  occur.  Because  the  impact-generated  plasma  has  an 
approximate  mean  energy  of  10-120  eV  only,  the  plasma  electron 
temperature  is  well  below  the  critical  temperature  T*  for  negative 
charging  to  occur.  Thus,  no  negative  charging  onset  would  occur 
for  the  spacecraft  interacting  with  the  plasma  generated  by  meteor 
impact.  Because  the  spacecraft  charging  theorem  is  independent  of 
the  plasma  density  n,  the  preceding  result  is  independent  of  how 
much  plasma  is  generated  on  meteor  impact. 

Sudden  and  Net  Charge  Emission 

If  the  charge  produced  on  neutral  particle  impact  is  not  neutral, 
meaning  there  is  a  net  charge  of  one  sign  produced  from  the  impact 
site  and  its  neighborhood  (Fig.  6),  charging  would  result.  Earlier 
measurements  of  separate  charge  production  have  been  reviewed 
by  Pailer  and  Griin.*''’  The  authors  concluded*'''  that  a  typical  value 
for  electron  yield  due  to  neutral  particle  impact  is  0.2-0.4,  rel¬ 
atively  independent  of  the  particle  energy,  provided  that  the  latter 
exceeds  the  ionization  threshold.  The  electron  (or  ion)  yield  is  the 
number  of  impact-generated  electrons  (or  ions)  emitted  from  or  near 
the  impact  site  per  impacting  particle.  Their  ion  yield*’*'  seems  in¬ 
conclusive,  varying  widely  from  0.01  to  2.  Later  measurements*'*’*'® 
confirmed  the  value  0.2  for  Yg  and  reported  0.04  for  T, ,  the  yield  of 
ions.  Schmidt  and  Arends*'®  used  aluminum,  gold,  and  while  paint  as 
targets  and  Xe,  Ar,  Kr,  and  air  molecules  as  beam  species,  whereas 
Rudenaucr  and  Steiger*'*  used  Xe  and  Ar  as  neutral  particle  beams 
and  an  aluminum  alloy  as  target.  Whereas  these  results  are  indica¬ 
tive  of  higher  eleep'on  yield  than  that  of  ions,  the  energies  of  the 
electrons  and  ions  have  not  been  measured.  The  yields  and  energies 
are  important  for  determining  the  level  of  charging,  as  discussed  in 
the  next  paragraph. 

Even  a  small  imbalance  of  electron  and  ion  yields  would  produce 
a  substantial  net  charge  on  impact  by  a  meteoroid,  if  the  yields  are 
large.  The  spacecraft  potential  is  determined  by  current  balance. 
At  geosynclu-onous  altitudes,  the  ambient  current  is  of  the  order 
of  milliampcre.'®’^*  For  an  emitted  electron  current  exceeding  the 
ambient  currents,  the  former  controls  the  current  balance,  driving 
the  spacecraft  positive  in  voltage.  If,  however,  the  emitted  electrons 
have  a  few  electron  volts  only,  the  charging  level  cannot  exceed  a 
few  volts. 

Our  conjecture  is  that  the  higher  electron  yield  (Yg  >  T,  )  is  due  to 
the  secondary  electrons  coming  out  from  beneath  the  surface  of  the 
target  material  in  the  neighborhood  of  the  impact  crater.  Secondary 
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Fig.  6  Plasma  and  neutral  gas  generated  by  the  hypervelocity  impact  of  a  particle  on  a  solid  surface. 
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electrons  have  typically  a  few  electron  volts  in  energy.  Charging 
cannot  exceed  a  few  volts  if  the  emitted  electrons  have  only  a  few 
electron  volts. 


We  now  calculate  the  quantity  (i.  Let  a  ionizations  be  generated 
by  an  electron  traveling  through  a  distance  dx.  For  n  electrons  at  x, 
the  number  of  ionizations  at  x  +  dx  is  given  by 


Triggering  Spacecraft  Discharges 

As  discussed  already,  hypcrvclocity  particulate  impact  on  space¬ 
craft  may  produce  1 )  mechanical  damage  to  the  spacecraft  and  2) 
plasma  and  neutral  vapor  that  pose  a  hazard  to  electrical  systems. 
These  effects  arc  not  mutually  exclusive  and  may  occur  at  the  same 
time;  next  is  a  brief  presentation  of  scenarios  where  these  hazards 
can  play  a  role. 

A  potential  weak  point  susceptible  to  hypcrvclocity  impact  haz¬ 
ards  arc  the  exposed  cables  on  a  spacecraft  (scenario  1).  The  cable 
sheath  may  have  weathered  as  a  result  of  years  of  operation  in  the 
harsh  space  environment.  Depending  on  the  energy  of  the  impacting 
particle  and  the  strength  of  the  target  material,  the  impacting  particle 
may  penetrate  or  even  tear  apart  the  worn  cable  sheath  and  cause 
electrical  disruption  to  the  circuits.  A  case  in  point  is  the  failure 
of  MSTI-2  satellite,  which  was  attributed  to  the  impact  of  orbital 
debris  on  a  wire  bundle.'*  ''^ 

Another  potential  weak  point  pertaining  to  this  type  of  problem 
is  the  solar  panel  power  system,'*®  as  illustrated  in  Fig.  7  (scenario 
2).  We  have  shown  that,  over  the  course  of  a  satellite  lifetime,  the 
probability  for  meteor  impacts  of  solar  panels  to  depths  of  5  mm  is 
quite  high.  Such  a  depth  is  sufficient  to  cause  great  damage  to  the 
solar  panels.  A  large  Leonids  meteoroid,  though  rare,  may  penetrate 
deeper.  A  case  in  point  is  the  failure  of  the  Olympus  satellite.  It 
occurred  suddenly  during  the  time  of  maximum  activity  of  a  Perseids 
shower.  The  lime  was  in  a  gcomagnetically  quiet  period,  which  rules 
out  surface  and  deep  dielectric  charging  as  a  cause.  The  gyro  and 
its  power  supply  failed.  The  cause  was  attributed  to  a  meteoroid  hit 
on  the  solar  panel  system  or  a  cable  connecting  the  system  to  the 
gyro.'®''”’'**’ 

If  the  impact  occurs  near  solar  cells  that  arc  at  different  potentials, 
the  impact-generated  plasma  might  then  connect  the  different  solar 
cells  (scenario  3).  The  number  of  electrons  is  maintained  constant  or 
increases  in  time,  provided  that  the  number  of  electrons  m  generated 
between  a  cathode  and  an  anode  sustains  or  increases  in  time.  This 
process  is  not  limited  to  solar  cells  (Fig.  7)  but  is  applicable  to  two 
or  more  nearby  differentially  charged  configurations.  We  discuss 
now  two  scenarios  in  which  the  number  of  electrons  is  maintained 
constant  or  increases  in  time. 

Let  /L  ionization  events  be  generated  in  the  neutral  gas  by  an 
electron  transit  between  a  cathode  and  an  anode.  Both  the  electrons 
and  the  neutral  gas  arc  produced  by  a  meteor  impact  on  a  surface. 
The  cathode  and  anode  can  be  two  nearby  differentially  charged 
surfaces.  The  ji  electrons  arc  accelerated  toward  the  anode,  and 
the  fi  ions  arc  accelerated  toward  the  cathode.  If  an  ion  impact  on 
the  cathode  can  produce  y  electrons,  then  n  ion  impacts  would 
produce  fiy  electrons  from  the  cathode.  A  necessary  condition  to 
sustain  the  discharge  is  the  inequality 

^ly  >  1  (12) 

The  newly  created  electron  would  start  its  journey  toward  the  anode, 
thereby  creating  at  least  /r  new  electron-ion  pairs  to  sustain  a  low- 
level  discharge  or  initiate  an  avalanche  ionization. 
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Fig.  7  Electrons,  ions,  and  neutral  gas  generated  by  a  hypervelocity 
impact  on  a  solar  panel. 


dn  =  na  dx  (13) 

The  total  number  of  ionizations  from  the  cathode  (x  =  0)  to  an  anode 
(x  =  d)  is  given  by 


n{d)  =  n(0)  expfarf) 


(14) 


where  /i(0)  is  the  initial  number  of  electrons  starting  the  journey  at 
X  =  0.  Thus,  fi  =  expiad). 

For  electron  impact  ionization  of  neutral  gas,  a  is  given  by^** 


a  <xN  j  dEEi  f{E)a{E)  (15) 

where  N  is  neutral  density  and  E  the  electron  energy.  Thus, 


=  exp  Nd  J 


dEEif(E)aiE) 


(16) 


Applying  the  result  [Eq.  (16)]  to  discharges  in  meteor-generated 
plasma,  we  sec  that  fi  is  proportional  to  cxp(Al)  and  depends  on  the 
integral  of  the  electron  density  distribution  function  /(£).  Thus, 
both  the  neutral  gas  density  and  the  plasma  electron  distribution 
function  play  important  roles  in  mediating  a  discharge.  The  effect 
of  a  sustained  low-level  discharge  can  drain  currents  and  degrade 
systems.  An  avalanche  discharge  can  short  a  circuit  or  cause  power 
supply  failure  rapidly. 

We  have  shown  the  roles  of  two  key  players,  namely,  the  plasma 
and  the  neutral  vapor,  both  generated  by  meteor  impacts  on  sur¬ 
faces.  The  inequality  [Eq.  (12)]  is  necessary  but  not  sufficient.  To 
be  sufficient,  one  needs  to  account  for  loss  mechanisms,  a  situation 
reminiscent  of  the  discharge  criteria  in  critical  ionization  velocity.^* 

A  heretofore  unsuspected  differentially  charged  configuration  is 
dielectric  surface  materials  impregnated  with  incoming  high-energy 
(mcgaelectron  volts  or  higher)  electrons  and  ions  (scenario  4).  Dur¬ 
ing  a  passage  of  a  high-energy  (megaelectron  volts  or  higher)  plasma 
cloud  or  a  solar  coronal  mass  ejection  (CME)  cloud,  the  energetic 
(mcgaelectron  volts  or  higher)  electrons  and  ions  penetrate  into 
materials  to  different  depths  and  stay  there,  if  the  conductivity  is 
low.  For  example,  megaelectron  volt  electrons  penetrate  to  10^  mil 
in  aluminum,  whereas  megaelectron  volt  protons  to  4  x  10"*  mil 
only.^  After  many  hours  or  days,  a  double  layer  is  formed  inside 
the  dielectric.  The  electron  flux  in  space  is  two  orders  of  magni¬ 
tude  greater  than  that  of  ions.'*^  '*®  However,  with  a  predominantly 
electron  layer  formed  inside,  the  dielectric  surface  and  its  vicin¬ 
ity  experience  an  electric  field  that  attracts  ions.  When  the  ambient 
environment  becomes  quiet  after  days  of  energetic  cloud  passage, 
the  ambient  plasma  returns  to  low  energy  again.  Low-energy  ions 
can  be  attracted  by  the  surface  electric  field,  thus  building  up  the 
shallow  ion  layer.  Nature  tends  to  evolve  toward  neutralization  and 
equilibrium.  Even  if  the  low-energy  ions  can  completely  neutralize 
the  surface  electric  field,  they  cannot  recombine  with  the  deeply 
layered  electrons  because  they  are  separated  at  a  distance.  Thus,  a 
nonequilibrium  configuration  ensues  and  is  unstable.  Shorting  the 
two  layers  would  lead  to  a  discharge,  which  is  called  an  anodized 
discharge.^^  A  hypervelocity  meteor  impact  penetration  can  trig¬ 
ger  such  shorting  (Fig.  8)  because  the  ionization  track  made  by  the 
meteor  impact  can  provide  a  conduction  channel  for  a  discharge 
between  the  electrons  and  ions  staying  at  different  depths.  A  case 
in  point  is  the  recent  failures  of  three  geosynchronous  communica¬ 
tion  satellites,  namely,  ANTK-1,  AT&T  Telstar-401,  and  Motorola 
Galaxy-4.  Although  the  true  causes  of  their  failures  may  never  be 
known,  all  three  of  them  occurred  after  (but  not  during)  the  passage 
of  a  CME. 

Another  hazardous  scenario  combining  the  penetration  and 
plasma/vapor  production  properties  of  hypervelocity  particulate  im¬ 
pact  is  as  follows  (scenario  5).  Depending  on  the  momentum  of  the 
meteoroid  or  particulate  and  the  thickness/strength  of  the  wall,  com¬ 
plete  penetration  may  occur.  If  the  penetration  goes  through  a  wall. 
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Fig.  8  Sequence  (a,  b,  c,  d)  of  double-layer  formation  in  a  dielectric 
material. 


Fig.  9  Meteoroid  penetrating  through  the  spacecraft  casing. 


the  plasma  and  neutral  vapor  can  come  out  from  the  other  side  of 
the  wall.  This  would  have  important  consequence  for  the  electron¬ 
ics  inside.  The  exposed  wires  inside  and  the  circuit  boards  may 
be  shorted.  The  alkali  vapor  released  from  the  meteoroid  may  re¬ 
act  chemically  with  certain  materials  such  as  aluminum  (Fig.  9). 
For  example,  sodium,  which  is  known  to  be  reactive,  can  react  with 
aluminum-bearing  minerals  to  produce  the  sodium  aluminosilicate^^ 
NaAlSisOg. 

The  generation  of  plasma  by  a  sudden  impact  can  produce  elec¬ 
tromagnetic  waves,^'*  which  then  can  introduce  noise  and  spurious 
signals  into  the  antenna  system^^^^*  (scenario  6).  Finally,  we  briefly 
point  out  that,  as  an  impact  occurs,  the  rapid  burst  of  neutral  gas,  elec¬ 
trons,  and  ions  would  give  off  a  flash  of  optical  and  electromagnetic 
radiation  signatures.  The  optical  lines  of  the  flash  would  depend  on 
the  compositions  of  the  spacecraft  surface  and  the  meteoroid  or  de¬ 
bris.  Collisions,  excitation,  ionization,  and  recombination  occur  in  a 
short  period  before  diffusion  transport  takes  over.  The  optical  emis¬ 
sion  and  even  the  impact-generated  dust  cloud  itself  can  temporarily 
degrade  the  performance  of  optical  instruments  onboard.  Dust  at¬ 
tachment  or  molecular  adsorption  on  optical  instruments  could  be 
hazardous.  A  more  detailed  discussion  of  this  subject  is  beyond  the 
scope  of  this  paper. 


Mitigation 

Meteors  are  rarely  taken  into  consideration  when  designing  space¬ 
craft.  It  is  simply  assumed  that  the  risk  is  low.  There  are,  however, 
a  few  mitigation  methods  that  have  been  employed  in  the  past.  The 
most  common  methods  for  reducing  the  impact  of  meteoroids  on 
spacecraft  are  preventive  in  nature.  These  methods  consist  of  1 )  min¬ 
imization  of  the  probability  of  meteoroid  impacts,  2)  minimization 
of  the  risk  of  a  short  circuit  or  a  discharge  should  an  impact  occur, 
and  3)  protection  of  sensitive  electronics  behind  some  shields  to 
avoid  impact  penetration. 
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Fig.  10  Whipple  shield. 


In  anticipation  of  an  intense  meteor  shower,  it  is  possible  to  orient 
the  solar  panels  in  a  direction  parallel  to  the  meteor  stream.  This 
minimizes  the  effective  area  normal  to  the  meteor  stream  velocity 
so  that  the  probability  of  a  strike  is  minimized.  It  is  also  possible  to 
turn  sensitive  instruments,  such  as  optics,  away  from  the  direction 
in  which  shower  meteors  will  come.  Another  common  practice  is 
to  shut  off  all  of  the  nonessential  electrical  power  on  a  spacecraft, 
during  the  period  of  an  intense  meteor  shower,  to  minimize  the 
probability  of  a  short  circuit  or  a  discharge  in  case  an  impact  occurs. 
Moving  a  spacecraft  to  the  side  of  the  Earth  opposite  the  shower 
radiant,  thereby  using  the  Earth  as  a  shield,  would  also  be  helpful, 
but  moving  a  spacecraft  is  often  a  difficult  task. 

Whipple^^  suggested  a  double- wall  shield,  comprising  a  bumper 
shield  and  a  primary  shield.  If  an  impinging  particle  penetrates 
through  the  bumper  shield,  the  molten  particle  or  vapor  generated, 
together  with  the  small  debris  from  the  bumper,  will  impact  the 
primary  shield  with  dispersed  instead  of  direct  energy  (Fig.  10). 
This  provides  a  better  opportunity  for  shielding.  Modem  designs 
of  multiwall  systems  using  high-strength  composite  materials  are 
improvements  over  the  original  Whipple  design  (see  Ref.  60).  Note, 
however,  that  few  satellite  employ  the  Whipple  shield. 

Conclusions 

We  have  emphasized  the  hazards  of  hypervelocity  particle  im¬ 
pacts  in  space.  Hypervelocity  particles  generate  great  energies  on 
impact  on  spacecraft  surfaces.  Based  on  other  laboratory  measure¬ 
ments  of  impact  penetration  depth,  empirical  formulas  of  penetra¬ 
tion,  and  plasma  production,  we  have  pointed  out  the  importance 
of  charge  imbalance  in  the  plasma  produced  on  impact.  As  shown 
in  the  discussion,  charge  imbalance  can  cause  spacecraft  charg¬ 
ing.  We  have  also  discussed  some  worst-case  scenarios  that  may 
arise  from  hypervelocity  impact  on  spacecraft.  To  sustain  a  diode- 
type  discharge,  it  is  necessary  to  satisfy  a  simple  Townsend-type 
inequality.  The  sufficient  conditions  would  require  accounts  of  loss 
mechanisms.  An  optical  or  electromagnetic  flash,  with  characteristic 
signature  lines,  may  result  from  hypervelocity  impact. 
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